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Abstract 

Moments of the hadronic invariant mass and of the lepton energy spectra 
in semileptonic B decays have been determined with the data recorded by 
the DELPHI detector at LEP. From measurements of the inclusive 6-hadron 
semileptonic decays, and imposing constraints from other measurements on b- 
and c-quark masses, the first three moments of the lepton energy distribution 
and of the hadronic mass distribution, have been used to determine parameters 
which enter into the extraction of |V c fe| from the measurement of the inclusive 
6-hadron semileptonic decay width. The values obtained in the kinetic scheme 
are: 



m 6 (l GeV) = 4.591 ± 0.062 ± 0.039 ± 0.005 GeV/c 2 

m c (l GeV) = 1.170 ± 0.093 ± 0.055 ± 0.005 GeV/c 2 

/4(1 GeV) = 0.399 ± 0.048 ± 0.034 ± 0.087 GeV 2 

p 3 D = 0.053 ±0.017 ±0.011 ±0.026 GeV 3 , 

and include corrections at order l/mf. 

Using these results, and present measurements of the inclusive semileptonic 
decay partial width of 6-hadrons at LEP, an accurate determination of |V c f>| is 
obtained: 

|V c6 | = 0.0421 x (1 ± 0.014 meas . ± 0.014 fit ± 0.015^.) . 
(Accepted by Euro. Phys. Journ.) 
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1 Introduction 

Several years ago it was proposed to obtain an accurate value of \V c b\ by comparing the 
measurement of the inclusive semileptonic decay partial width in the process b — > ci~V£, 
with the corresponding theoretical expression, obtained using the Operator Product Ex- 
pansion (OPE) formalism, applied in the heavy quark mass limit [1]. A recent appraisal 
of the limitations of this approach can be found in [2] from which the following expression, 
for the semileptonic decay width, has been taken: 

T sl (b^ci-V e ) = \V cb \ 2 (l + A ew )A^(r,p) 



z (r 



2m 2 b (p) 



2(1 - r f M 2 + d(r)^4 + 

mftfj,) m 3 b (p) 



(1) 



In this expression, zo(r) is the tree-level phase space factor and r = m 2 (p) / 'm 2 (p) . Def- 
initions for the other quantities are given in [2]. Equation (1) contains an expansion 
in a s , corresponding to perturbative QCD corrections expressed in A pert (r,p), and an 
expansion in l/m&, corresponding to non-perturbative QCD contributions. An auxiliary 
scale /i(= 1 GeV) is introduced to demark the border between long- and short-distance 
dynamics in OPE. Numerically, A ew , corresponding to the ultraviolet renormalization of 
the Fermi interaction, is well-known and amounts to +1.4% [3]; A pert (r,p) corresponds 
to ~ —9% corrections and the non-perturbative terms have typically few percent contri- 
butions [2]. The smallness of these last corrections comes, partly, from the fact that the 
term proportional to l/m& is absent [4,5] in Equation (1). The quantities p 2 , Pq, p% and 
p\ s denote the expectation values of the kinetic, chromomagnetic, Darwin and spin-orbit 
operators respectively. These parameters have to be determined by experiment. From the 
mass splitting between B* and B meson the following value Pq(1 GeV) = (0.35loQ2)GeV 2 
has been obtained [6]. The value for p 2 is less certain; in this regularization scheme the 
inequality p 2 (p) > Pg(p) holds for any normalization scale. Constraints have been also 
established for terms contributing at order 1/mjj. The Darwin parameter, p%(p), must be 
positive and the spin-orbit operator value, p|,s(/x), * s expected to be negative and to sat- 
isfy —p\s — 9% [2]- I n addition, it has been demonstrated that the value of T s i(b — > cl~Vg) 
is rather insensitive to the exact value of p\ s [7]. 

A few lessons can be drawn from the previous considerations: 

• the largest correction, to the naive free quark decay model, is expected to originate 
from perturbative QCD. Its evaluation is closely connected to the definition adopted 
for the running quark masses; 

• in addition to the values of heavy quark masses, m& and m c , two parameters need 
to be determined by experiment, p 2 , and p 3 D , to have a control of non-perturbative 
QCD corrections up to the l/mf order; 

• the definition for the quark masses has to be consistent with both perturbative and 
non-perturbative dynamics. 

It is then expected that the value of \V c b\, extracted in this way from inclusive b- 
hadron semileptonic decays, can be determined with a relative uncertainty from theory 
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at the 1.5% level (see [2] for a detailed breakdown of contributing sources in this num- 
ber). To match this accuracy, an experimental control of the parameters governing non- 
perturbative QCD corrections, even at a modest level, is required and is the main purpose 
of the present analysis. 

The determination of moments of the lepton energy and of the hadronic mass spectra 
in B — * X c £ve decays provides important information on these parameters since they can 
be analysed using the same formalism, and since corresponding expressions depend on 
the same non-perturbative parameters entering in Equation (1). The following notations 
for moments of the lepton energy spectrum have been used: 

M[ =< E\ > and =< (E* e - < E\ >) n >, n > 1, (2) 

where E\ denotes the value of the lepton energy in the 6-hadron rest frame; and for 
moments of the hadronic mass system: 

K =< (m 2 H - m% m Y > and Mf =< (m\- < m\ >)« > . (3) 

fn sp i n = 1.97375 GeV/c 2 denotes the spin averaged D meson mass, which is equal to the 
weighted average of the D and D* masses. The OPE expresses lepton moments through 
quark masses as a double expansion in a s and l/m^. 



2 2 3 3 

7i rn h ™i ml ml 



The higher coefficient functions b n (r), c n (r), ... are also perturbative series in a s . Due to 
the kinematic definition of the hadronic invariant mass, Mx, the general expression for 
the hadronic moments includes explicitly the value of the 6-hadron mass, M#: 



1=0 



M B - m b 
m b 



\E nl (r) + a nl (r)^ + b nl (r)^ + c nl (r)^ 



n mi mi 

+ dm(r)4 + *m(r)^f + ...V (4) 
ml ml J 

Numerical values for all r-dependent functions entering into these expressions can be 
found in [7]. 

In the following, the //-scale independent third order correction term, p 3 D has been 
fitted in place of p 3 D . The two quantities are related through the expression: p 3 D pa 
p 3 D (l GeV) -0.1 GeV 3 [8]. 

The consistency of the bounds set by moments of different distributions, and with 
other data, tests the underlying theory assumptions. 

It can be noted that the approach used here does not rely on the validity of an expan- 
sion in l/m c , as already advocated in [7]. The value of the charm quark mass, entering 
into the expression of moments, is taken as a parameter whose value has also been fitted 
using data. 

Measurements of moments have been reported by the CLEO [9-11], BaBar [12-14] 
and BELLE [15,16] Collaborations operating at the T(4S) resonance, and by DELPHI 
using preliminary data taken at the Z pole energy [7]. Results have also been recently 
published by CDF [17]. 

While there is an obvious advantage in measuring the energy spectra in events where 
the decaying B rest frame almost coincides with the laboratory frame, low energy particles 
cannot be identified there. It is thus necessary to rely on models for extrapolating the 
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lepton energy spectrum to zero energy or to resort to computations for a truncated 
spectrum. Performing this analysis at energies around the Z peak ensures sensitivity 
to almost the full lepton spectrum, thus reducing modelling assumptions. The main 
challenge at the higher energy is, for the lepton energy moments analysis, the accurate 
determination of the B rest frame. 

This paper presents the results obtained from analyses of the data recorded with 
the DELPHI detector at LEP on moments of the hadronic mass and charged lepton 
energy distributions. The analysis procedures are discussed in Sections 4 and 5. In 
Section 6 these results are then used as inputs of a multi-parameter fit to determine 
the relevant corrections contributing at 0(l/mfy and 0(1/ m^), together with the heavy 
quark masses. The use of higher moments guarantees a sensitivity to these parameters 
and the simultaneous use of the hadronic and leptonic spectra ensures that a larger 
number of parameters can be kept free in the fit. We discuss the results both in terms 
of the extraction of the parameters and the implications for jV^,], and as a consistency 
check of the underlying theoretical assumptions. 

In addition, production and decay properties of broad D** states have been studied. 
They are reported in Section 3. D** refers, in the present analysis, to all hadronic 
systems of mass higher than the D*. 6-meson semileptonic decays with charmed hadrons 
emitted in the final state correspond to B — > D, D*, D**£~z7^ transitions. D** states 
are about 30% of 6-hadron semileptonic decays [18,19]. They can be resonant or non- 
resonant hadronic D^rar systems. Resonant states are supposed to be dominated by 
L=l, orbitally excited cq states. There are four such states: Dq, D*, DJj and D± with 
respectively J p = + , 1 + , 2 + and 1 + . The two 1 + states, having common final states, 
mix and physical states are expected to decay into almost pure D or S wave D*7r or 
Dp final states [20]. The two states (D£ and D^ which, because of angular momentum 
and parity conservation, have to decay into a D wave are expected to be narrow and 
their measured widths are of the order of 20 MeV/c 2 [21]. The contribut ion from these 
narrow resonances has been measured [19,22] and it amounts to about one third of all 
D** states. Non-resonant or broad resonant states are thus expected to have a dominant 
contribution in 6-hadron semileptonic decays. In addition to L=l mesons, there could 
be states corresponding to higher angular momentum values or to radial excitations. All 
such states are expected to be broad [23]. Contributions from non-resonant Dtt or D*7r 
final states have been evaluated in [24] and this model has been considered in the present 
analysis. 

In a previous DELPHI publication [25], production fractions of D** states in 6-hadron 
semileptonic decays have been measured. The main focus of the present analysis is on the 
mass distribution of these states. This was made possible by improving the separation 
between signal and background events with respect to the previous analysis. 

2 Data analysis 

This study is based on 6-hadron semileptonic decays recorded with the DELPHI detec- 
tor at LEP from 1992 to 1995. Since the determination of the moments of the hadronic 
mass and charged lepton energy distributions have different requirements, two analyses 
have been performed. The first focuses on the exclusive reconstruction of D** states. 
For the second analysis a lepton sample with low background is required. This section 
presents those parts of the data selection and event reconstruction procedure which are 
common to both analyses. 
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2.1 Hadronic event selection and simulation 

In order to select hadronic Z decays, standard hadronic selection cuts have been ap- 
plied. Each event has been divided into two opposite hemispheres by a plane orthogonal 
to the thrust axis. The polar angle 1 of the thrust axis of the event had to satisfy the 
requirement |cos#| < 0.95. Charged and neutral particles have been clustered into jets 
using the LUCLUS [26] algorithm with the resolution parameter dj i n = 5 GeV/ c. About 
3.4 million events have been selected from the LEP1 data sets. 

The JETSET 7.3 Parton Shower [26] program has been used to generate hadronic Z 
decays, which were passed through the detailed detector simulation DELSIM [27] and 
processed by the same analysis chain as the data. A sample of about nine million Z — > qq 
events has been used. To increase the simulation statistics, an additional sample of about 
3.6 million Z — > bb events, equivalent to about 17 million hadronic Z decays, has also been 
used. Statistics of the analysed hadronic samples are given in Table 1. 



Year 


Real data 


Simulated 
Z^qq 


Simulated 
Z -> bb 


1992+1993 
1994+1995 


1355805 
2012921 


3916050 
5012881 


1096199 
2495335 


Total 


3368726 


8928931 


3591534 



Table 1: Analysed number of events. In 1992 and 1993 only two-dimensional vertex 
reconstruction was available. 

Z — > bb events have been selected using an event 6-tagging technique [28] based on the 
reconstructed impact parameters of particle tracks. 

Events for the exclusive analysis have been selected by requiring the presence of one 
tagged lepton candidate with momentum p > 2 GeV/c and of a D°, D + or D* + candidate 2 
in the same event hemisphere. 

In order to measure moments of the lepton energy distribution in inclusive b — > X c lVi 
decays, events have been required to contain one tagged lepton candidate with momentum 
p > 2.5 GeV/c (for muons) or p > 3 GeV/c (for electrons). 

2.2 Muon identification 

Muons have been identified based on the response of the Muon Chambers. Details can 
be found in [27]. 

For the inclusive lepton analysis muon candidates have been accepted if they fulfilled 
the "standard" selection criteria, their momenta in the lab frame exceeded 2.5 GeV/c and 
were contained within the polar angle intervals: | cos# M | < 0.62 or 0.68 < | cos# M | < 0.94, 
defining the barrel and the forward regions. The muon identification efficiency has been 
measured with Z — > events, in the decays r — > fiu T V ^ and in two-photon 77 — ■> fi + fi" 
events. A mean efficiency within the acceptance region of 0.82 + 0.01 has been found, with 
little dependence on the muon momentum and on the track polar angle. This agrees with 
simulation, both in absolute value and in the momentum dependence, within a precision 
of 2%. 

1 In the DELPHI coordinate system, z is along the electron beam direction, <j> and R are the azimuthal angle and radius 
in the xy plane, and 8 is the polar angle with respect to the z axis. 

Throughout this paper charge-conjugate states are implicitly included. 
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The probability for a hadron to fake a muon has been estimated on anti-6 tagged events. 
After subtracting the expected remaining muon content in this sample, the misidentifi- 
cation probability for hadrons has been found to be (0.52 ± 0.03)% in the barrel and 
(0.36 ± 0.06)% in the forward regions, respectively. Applying the same procedure to sim- 
ulation events gave however values lower by factors of 2.03 ± 0.12 in the barrel and of 
1.22 ± 0.20 in the forward regions respectively. The simulation predictions have therefore 
been corrected for these factors. 

In the exclusive D** analysis, muon candidates have been accepted if they fulfilled 
the "loose" selection criteria and their momenta exceeded 2 GeV/c. The corresponding 
efficiency is ~ 80% and the hadron misidentification probability is ~ 1%. With this 
selection, the correction factors to be applied to simulated events in which the candidate 
lepton is a misidentified hadron have been found to be 1.44 ±0.05 (1.61 ±0.05) in 1992-93 
(1994-95) data samples 

2.3 Electron identification 

Electron candidates have been tagged within the range 0.03 < |cos# e | < 0.72 in 
polar angle, based on the combination of the response of the HPC, the specific ionization 
dE/dx in the TPC and the RICH Cherenkov detector [27]. A momentum dependent cut 
applied on the neural net output variable, which provides a constant efficiency over the 
full momentum range, has been applied. 

For the inclusive lepton analysis, electron candidates have been selected with a selec- 
tion cut corresponding to 65% efficiency and requiring momenta greater than 3 GeV/c. 
The probability for a hadron to fake an electron is about 0.4%. The electron identifica- 
tion efficiency has been measured from data by means of a sample of isolated electrons 
extracted from selected Compton events and one of electrons originating from photon 
conversions in the detector. The ratio between the efficiencies measured in data and 
simulated events has been parametrized as a function of the transverse momentum and 
polar angle of the particle track. Results are summarised in Table 2. A corresponding 
correction factor has been applied to simulated qq events. 

The probability for mis-tagging a hadron as an electron has also been measured using 
data, by selecting an anti 6-tagged background sample, as for the muons. 

Electrons from photon conversions, mainly produced in the outer ID wall and in the 
inner TPC frame, have been rejected by removing electron candidates originating at a 
secondary vertex and carrying little transverse momentum relative to the direction defined 
from the primary to the secondary vertex. The ratio of the measured misidentification 
probability in data to that in simulated events is given in Table 2. 

In the exclusive D** analysis, electron candidates have been selected with 75% effi- 
ciency and requiring momenta greater than 2 GeV/c. The probability for a hadron to 
fake an electron is about 1%. Correction factors applied to simulated events are also 
given in Table 2. 

2.4 Hadronic decay reconstruction 

The reconstruction of D° and D* + candidates, in which the D° decays into K~7r + , 
K~7r + 7r + 7r~ or K~7r + (7r°), is explained in detail in [29]. The reconstruction of the D + — > 
K~7T + 7r + decay is based on a similar approach. For all decay channels, the main steps of 
the analysis consist of: 

• reconstructing a D decay vertex from its charged decay products; 
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1992 


1993 


1994 


1995 


65% efficiency cut 
Efficiency (data/MC) 


0.83±0.02 


0.83±0.02 


0.92±0.02 


0.93±0.02 


Misid. Prob. (data/MC) 


0.57±0.04 


0.77±0.05 


0.76±0.05 


0.70±0.06 


75% efficiency cut 
Efficiency (data/MC) 


0.89±0.02 


0.88±0.02 


0.94±0.02 


0.94±0.02 


Misid. Prob. (data/MC) 


0.61±0.04 


0.77±0.03 


0.80±0.02 


0.74±0.03 



Table 2: The ratio between values measured in real and simulated events for electron 
identification efficiency and probability of tagging a hadron as an electron 



• selecting a mass window centred on the nominal D mass; 

• reconstructing a B decay vertex using the D trajectory and a charged lepton; 

• requiring a minimum distance between the B decay vertex and the main vertex and 
also between the B and the D decay vertices; 

• imposing a minimum momentum on the D (6 GeV/ c) and on the D-£ (10 GeV/ c) 
candidates; 

• requiring a D-£ mass between 2.5 and 5.5 GeV/c 2 . 

3 D** production in 6-hadron semileptonic decays 

In the following sections the exclusive analysis leading to the measurement of hadronic 
moments is presented. D°, D + , D* + and D** reconstruction is described in Section 3.1. 
In Section 3.2 a discriminant variable that has different sensitivity to D** signal and 
background events is defined. This variable is used in Section 3.3 to measure the amount 
of D**£~Ve states in the data. From the study of the D** mass distribution, branching 
ratios and properties of different D** states are measured in Sections 3.4 and 3.5. 

3.1 D°, D+, D*+ and D** reconstruction 

Mass distributions of D° and D + candidates and distributions of the mass difference 
between the D°7r + and the D° candidate, in the case of channels involving a D* + , have 
been used to define signal and sidebands samples (Figure 1). Events from the sidebands 
have been used to evaluate the level of the combinatorial background under the D or D* + 
signal. 

For D** states, decay channels into a D'*' and, at most, two pions have been considered. 
When searching for D** — > D°7r + decays, a veto has been applied against the D* + by 
removing candidates with a mass difference m(D°7r + ) — m(D°) < 0.1465 GeV/c 2 , for D° 
mesons decaying into charged particles only 3 . As the combinatorial background under 
D signals is higher than for channels involving a D* + , values of the cuts given in [29] have 
been made tighter for D** — > Dn events, as compared with those applied to D** — > D* + 7r 
decays. 

Events have been selected by reconstructing a lepton, a and a charged pion whose 
trajectories are compatible with the hypothesis that they originate from a common sec- 
ondary vertex. D** decays considered in the present analysis involve always at least one 

3 For D°7r final states, only D° decaying into K~7r+ or K~7r+7r+7r~ have been used. 



7 




Figure 1: D° ; D + and D* + signals used in the present analysis, in which the D° meson 
decays into charged particles only and which correspond to events registered in 1992-95. 
Intervals used to define the signal and sideband regions are indicated. 
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charged particle track 4 emitted at the 6-decay vertex, in addition to the exclusively recon- 
structed D**) meson. The %** momentum has been required to be higher than 0.5 GeV/ c 
and the track must be associated to at least one measurement in the Vertex Detector. 

The overall efficiencies for selecting signal events (see Table 3), not including the decay 
branching fractions of D* + , D° or D + mesons into their considered decay channels, have 
been estimated from simulated events. Efficiencies are rather similar for the 92-93 and 
94-95 samples when considering channels with a D* + or with D° — > K~7r + in spite of the 
reduced performances of the VD for the 92-93 period. This is because, as the background 
level is rather low for these channels, loose cuts on the vertex separation have been 
applied. For the other channels, which require tighter cuts, efficiencies are markedly 
lower in 92-93. These values have already been corrected for differences between the 
actual measured lifetimes of 6-hadrons and that used in the simulation. 



decay channel 


92-93 MC 


94-95 MC 


D** -> D U 7T + 






D u -> K-7T+ 

D° -> K-TT+TT+TT- 


(10.4 ±0.5)% 
(3.0 ±0.3)% 


(13.0 ±0.5)% 
(4.8 ±0.3)% 


D** -> D+7T- 






D+ -> K-7T + 7T + 


(5.9 ±0.3)% 


(9.1 ±0.3)% 


D** -> D*+7T- 






D u -> K-7T+ 

D° -> K-TT + TT + TT- 
D° -> K-7T + (7T ) 


(12.0 ±0.3)% 
(5.0 ±0.2)% 
(6.6 ±0.6)% 


(13.9 ±0.4)% 
(5.8 ±0.2)% 
(7.3 ±0.6)% 



Table 3: Global efficiencies of the analysis chain to reconstruct and select simulated signal 
events. Quoted uncertainties are only statistical. In addition to the simulated events 
mentioned in Table 1, dedicated event samples corresponding to the different channels 
have been used to increase the statistics. 



3.2 Signal separation from background sources 

When considering D < -*' ) tc ± combinations, the main sources of background can be divided 
into two categories depending on whether they correspond or not to a real reconstructed 
D**) meson. The latter is the combinatorial background situated under the charm mass 
signal. Background events with a real D'*' can originate from the following sources: 

• the 7i~** candidate is not produced at the 6-decay vertex but comes from the beam 
interaction point (primary pion background); 

• the lepton originates from the weak decay of another charm particle emitted in the 
6-decay (cascade background); 

• the lepton originates from a r decay (tau background); 

• the candidate lepton is a misidentified hadron or a converted photon (fake lepton 
background) ; 

• the reconstructed charm meson originates from a cc event (charm background). 

A variable, used to isolate the signal from these backgrounds, has been defined from 
the probability distributions of several discriminant observables, whose shapes have been 

4 This particle is called it** in the following analysis . 
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obtained from the simulation and also directly from data, as in the case of the combina- 
torial background. The following observables have been used: 

• the lifetime-signed impact parameters of the 7r** relative to the main vertex of the 
event, in R0 and z projections, normalised to their uncertainty; 

• the normalised and lifetime-signed impact parameters of the n** relative to the 
secondary vertex, in R0 and z projections; 

• the normalised and lifetime-signed decay distance between the primary and the 
secondary vertices; 

• the cosine of the decay angle defined as the angle of the it** direction, boosted 
to the 6-hadron rest frame, relative to the B direction. This variable is uniformly 
distributed for the signal whereas it is peaked at negative values for backgrounds; 

• the x 2 probability for the secondary vertex, which should be uniformly distributed 
for the signal and peaked at small values for the main sources of background; 

• the two variables, d±, which depend on the presence of additional charged particles 
at the secondary vertex. They are defined in the following way: 

— all charged particles, other than the D** decay products and the lepton, emitted 
in the same event hemisphere as the 6-candidate, with momentum larger than 
500 MeV/ c, which form a mass with the D^tt**^ system lower than 6 GeV/ c 2 
and which have values for their impact parameters relative to the 6-decay vertex 
smaller than 2 and 1.5 a in R0 and z respectively, are considered; 

— selected particles, having the same (+) or the opposite (— ) charge as the lepton 
are considered separately. If there are several candidates in a class, the one with 
the largest impact parameter to the main vertex is retained and the quantity: 



is evaluated, where e and nsig are, respectively, the sign and the number of 
standard deviations of the track impact parameter relative to the main vertex 5 . 
The sign of the impact parameter is taken to be positive (negative) if the corre- 
sponding track trajectory intercepts the line of the jet axis from the main vertex 
downstream (upstream) from that vertex. 

As the track impact parameters can extend to very large values because of the 
relatively long decay time of 6-hadrons, the variables d± are taken to be equal to 
the logarithm of (1 + x±) and their sign is taken to be the same as x±. For events 
with no spectator track candidate, that is, with no additional tracks compatible 
with the 6-decay vertex, a fixed value of -4.0 is used for d±. For D^7r signal events, 
it is expected that no additional track is present at the 6-decay vertex whereas for 
]3(*) 7r + 7r - candidates, another track with a precisely defined charge correlation with 
the lepton, is expected. These properties have been used in the analysis in the 
following way. Selected events have been distributed into right-sign and wrong-sign 
candidates. Right-sign events correspond to D°7r + , D + 7r~ and D* + 7r~ pairs whereas 
wrong-sign events have an opposite sign pion. Since only D** + or D**° states can be 
produced in the semileptonic decay of a 6-hadron, where they are accompanied by a 
negatively charged lepton, wrong-sign combinations can receive contributions only 
from D(*)7r + 7r - final states. For right-sign combinations it has been required that 
signal events behave as if there is no additional charged particle track at the 6-decay 
vertex. This implies that, in the case of D^nn decays, only D°7r + 7r° decays or those 
involving two charged pions and where the n~ escapes detection, can contribute. For 

5 The second term of this equation is not considered for 1992-1993 data since the z coordinate is not measured. 




(5) 
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wrong-sign combinations, it has been required that signal events behave as if there 
was another charged particle track, of sign opposite to the n** and of trajectory 
compatible with the 6-decay vertex position. 

Probability distributions for each of these nine (seven for 1992-1993 data as the vertex 
detector measured only R0 coordinates) variables, have been obtained using the simula- 
tion for signal and background events. The agreement in the shape of the distributions of 
these variables for real and simulated events is illustrated in Figure 2-left in which the x 2 
probability distributions for secondary vertices for events selected in the sidebands of the 
signal have been compared. This distribution was selected as it is sensitive to possible 
differences between real and simulated events. 

The probabilities P S i gna i and Pbackg. have been obtained by multiplying the probability 
of each discriminating variable and a global discriminant R has been defined as: 

D P signal Pbackg. i r \ 

p Tp ' w 

1 signal < 1 backg. 

in which events corresponding to the two largest background sources, namely the com- 
binatorial and the primary pion background components, have been used to determine 
the probabilities for background events. The distribution of the variable R is peaked at 
+ 1(-1) for signal (background) events (see Figure 3). 

In Figure 2-right the distributions of this quantity for combinatorial background events 
selected in the sidebands of the D° signal for real and simulated events are compared. 
These distributions are in agreement but, to account for possible differences between real 
and simulated events and as the combinatorial background component is the largest of 
the background components for D°7r and D + 7r final states, real events situated in the 
sidebands of the signal have been used, in the following, to determine also the shape of 
the R distribution for the combinatorial background component. 



3.3 D** production rates 

Events corresponding to the D° — > K~7r + (7T°) decay channel have not been included 
in this study, as they are affected by a larger combinatorial background. 

The R variable distributions, for the different decay channels, have been measured 
separately for right- and wrong-sign samples and for the two data-taking periods. A 
binned likelihood fit [30] has been performed 6 in order to measure the branching fraction 
in each channel. The fractions of signal and primary pion background events have been 
left free to vary in the fit whereas the fractions of the other components with a real D**) 
have been fixed to values taken from external measurements. The discriminant variable 
distribution for combinatorial background events has been obtained from data using the 
sideband events. Distributions for the other background components have been taken 
from the simulation. 

Fractions which have been kept fixed in the fit have then been varied, in turn, according 
to their expected overall accuracy and the fit repeated to estimate the corresponding 
systematic errors. The fitting procedure has been verified on simulated events (see the 
last column of Tables 5-6). 

The simulation of cascade events comprises b — > D~DX events only. Corrections have 
therefore been applied to the rate and the topology of these events to account for the 
missing decay modes. Present measurements from ALEPH [31] and BaBar [32] on double 

6 Statistical uncertainties due to the finite number of analyzed simulated events are included in this fit. 
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Figure 2: Left: Distributions of the decimal logarithm of the \ 2 probability for secondary 
D°ir£~ vertices in events selected in the sidebands of the signal obtained using real (points 
with error bars) and simulated (histogram) events. Right: Distributions of the values of 
the discriminant variable corresponding to events selected in the sidebands of the signal 
for real data (points with error bars) and simulated events (shadowed histogram). The 
other histogram (hatched) corresponds to simulated events from combinatorial background 
situated in the signal mass region. The histograms have been normalized to the same 
number of entries. 
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Figure 3: Discriminant variable distributions for a) D* + 7r~ and b) D* + 7r + candidates. 
Points with error bars correspond to real events whereas the histograms show the fitted 
components. The fake lepton background in these plot includes r events. 

charm decays of 6-hadrons have been used, together with results on c-hadron inclusive 
and exclusive semileptonic decays from [21], to determine these corrections. They depend 
on the topology of the studied channel. For each channel, the correction factor on the 
branching fraction and the probabilities to have no track of same (V(0) ss ) or opposite-sign 
(V(0) os ) as the lepton have been evaluated (see Table 4). 



Channel 


Expected rate (%) 


corr. factor 


P(0)„(%) 


V(0)os(%) 


Right sign candidates 










b D u TT+rX, veto on D*+ 


1.61 ±0.29 


3.5 


33 ±4 


33 ±4 


b^D + n-£-X 


0.54 ±0.13 


2.3 


63 ±4 


15 ±3 


b -> D* + n-rX 


0.65 ±0.26 


2.7 


69 ±5 


13 ±3 


Wrong sign candidates 










b -> D u ir-£-X, veto on D*+ 


0.58 ±0.12 


2.6 


~ 1 


~ 


b -> D+n+rX 


0.43 ±0.09 


3.4 


13 ±2 


90 ±3 


b -» r>* + n + rx 


0.58 ±0.17 


4.2 


22 ±3 


89 ±2 



Table 4: Double-charm expected rates and probabilities for having no additional charged 
track in the considered decay channels. 

Measured branching fractions are summarized in Tables 5 and 6 respectively for right- 
and wrong-sign candidates. 

The following contributions to the systematic uncertainties have been considered: 

• uncertainties related to the values of external parameters, such as Rb, and the differ- 
ent branching fractions of charmed hadrons into the reconstructed final states. The 
values used in the present analysis, taken from [21], are given in Table 7; 
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Right sign candidates 


number of 


branching 


number of events 


Channel 


fitted events 


fraction (%) 


fitted (expected) MC 


b -> r> u TT + £-X, veto on D*+ 








D u -> K~-k + (92-93) 


109.0 ±27.5 


2.34 ± 0.59 


515.2 ± 67.7 (435) 


D° -> K-7T+ (94-95) 


60.3 ±20.3 


0.69 ±0.23 


939.6 ±67.5 (1001) 


D° -> K-7r + 7r + 7r- (92-93) 


13.2 ±20.0 


0.51 ±0.77 


262.8 ± 50.9 (243) 


D° -> K"7r + 7r + 7r" (94-95) 


59.0 ± 20.7 


0.93 ± 0.33 


774.5 ± 64.5 (756) 


Average 




0.89 ±0.18 












D+ -> K~n + n + (92-93) 


47.4 ± 18.8 


0.76 ±0.30 


350.1 ±47.6 (394) 


D+ -f K~n+n+ (94-95) 


36.2 ± 16.5 


0.25 ±0.11 


1101.1 ±63.0 (1132) 


Average 




0.31 ±0.10 












D u -> KTtt + (92-93) 


13.2 ± 9.3 


0.36 ± 0.26 


125.3 ±23.6 (128) 


D° -> KT7T+ (94-95) 


32.8 ± 8.8 


0.52 ±0.14 


282.6 ± 28.7 (303) 


D° -> K-tt+tt+tt- (92-93) 


16.4 ± 7.1 


0.55 ±0.24 


111.2 ± 20.2 (96) 


D° -> K-tt+tt+tt- (94-95) 


12.9 ± 7.2 


0.25 ±0.14 


239.5 ± 28.0 (227) 


Average 




0.40 ±0.09 





Table 5: Measured branching fractions for b — > Dtt£~X in right-sign combinations using 
the D decay branching fractions taken from Table 7. The last column gives the results 
obtained on simulated events. Numbers within parentheses correspond to the real number 
of simulated signal. 



• detector-dependent uncertainties such as those on the tracking efficiency, on the 
lepton identification efficiency and on the correction of differences between real and 
simulated events relative to the fake lepton rate; 

• differences between real and simulated events on track reconstruction accuracy. 
Distributions of the discriminant variable for signal-like events, obtained in real 
and simulated events, have also been compared using D* + — > D°ir + decays, with 
D° — > K~7r + or K~7r + 7r + 7r~. These events have been analysed using the same cri- 
teria as those applied to D** — > D°7r + candidates. A discriminant variable (R*) is 
constructed using the same input quantities as for the R variable, apart from d± 
whose effects have been evaluated separately. Possible differences between real and 
simulated events affecting track offset measurements, decay length reconstruction 
and x 2 probability of secondary vertices, in the case of a signal, can then be studied. 
Distributions of the R* variable obtained using 461 reconstructed D* + in real data 
and 4100 in simulated events are compared in Figure 4. The two distributions agree 
within the statistical uncertainties. A possible difference in the shape of these dis- 
tributions has been parametrized assuming a linear variation with the value of R*. 
The fitted slope is equal to —0.01 ± 0.08. In the following, the effect of a variation 
on the shape of the discriminant variable distribution, induced by a linear correction 
of slope equal to ±0.1, has been evaluated; 

• corrections of the cascade decay rate applied to simulated events have been evaluated 
with a 30% relative uncertainty (see Table 4). This variation is larger than most of 
the uncertainties attached to the correction applied to each channel and covers the 
dispersion of the values of these corrections. Variations in the cascade rates have 
been applied in a correlated manner to all channels when this was relevant. Effects 
from possible differences between real and simulated events for the d± distributions 
have been evaluated in [29] by studying the decay channel b — > D* + X£~V£. The 
probability to have no additional track, for signal events in which there is in reality 
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Wrong sign candidates 
Channel 


number of 
fitted events 


branching 
fraction (%) 


number of events 
fitted (expected) MC 


b -> D"tt-£-X, veto on D*+ 








D u -> K~tt+ (92-93) 
D° -> K--K+ (94-95) 
D° -> K-TT+Tr+Tr- (92-93) 

D° -» K-TT + TT + TT" (94-95) 


34.3 ± 17.8 
26.3 ± 15.3 
3.5 ±14.1 
1.4 ± 15.5 


0.74 ±0.38 
0.30 ±0.17 
0.14 ±0.56 
0.02 ±0.23 


57.0 ± 41.9 (70) 
136.8 ±36.6 (141) 

25.1 ±33.1 (40) 
143.1 ±41.1 (166) 


Average 




0.26 ±0.13 




6 -» D+TT+rX 








D+ K — TT+TT+ (92-93) 
D+ K — TT+TT+ (94-95) 


10.1 ± 17.2 

-9.9 ±15.8 


±0.16 ±0.27 
-0.07 ±0.11 


119.8 ±20.0 (87) 
190.6 ± 44.2 (226) 


Average 




-0.04 ±0.10 




b D*+TT+f— A 








D u -> K~tt+ (92-93) 
D° -> K~tt+ (94-95) 

D° -> K-TT + TT + TT" (92-93) 
D° -> K-TT + TT + TT" (94-95) 


-0.3 ±8.3 
0.62 ±6.3 
9.0 ±6.5 
-8.2 ±7.8 


-0.01 ± 0.22 
±0.01 ±0.10 
±0.30 ± 0.22 
-0.16 ±0.15 


not simul. 


Average 




0.00 ±0.07 





Table 6: Measured branching fractions for the different wrong-sign combinations. The 
last column gives the results obtained on simulated events and numbers within parentheses 
correspond to the real simulated signal. 



no other track coming from the 6-vertex, is of the order of 80% and we have taken a 
5% relative error. For events, in which there is at least one such track produced at 
the secondary vertex, the probability to miss it is estimated to be (20 ± 10%). For 
events selected in the 5m region, situated above the D* + signal, in real and simulated 
events, the probabilities for having no spectator track are of the order of 34% and 
differ by (2.5 ± 1.5)%. 

The uncertainty related to the finite statistics of simulated events has been included in 
the statistical uncertainty of fitted signal event numbers. 



parameter 


central value 


or hypothesis 


and uncertainty 


R* 


0.21664 ±0.00068 


BR(D*+ -> D°7T+) 


0.677 ±0.005 


BR(D° -> K-?T + ) 


0.0380 ±0.0009 


BR(D° -> K-7T + 7T + 7T-) 


0.0746 ±0.0031 


BR(D° -> K-7T + 7T ) 


0.131 ±0.009 


BR(D° -> K-£+u e ) 


0.070 ±0.003 


BR(D° -> K"K + ) 


0.00412 ±0.00014 


BR(D+ -> K-7T + 7T + ) 


0.091 ±0.006 




0.388 ±0.013 


P(b - BO) 


0.106 ±0.013 




1.542 ±0.016 ps 



Table 7: Values for the external parameters used in the analysis [21]. 
A summary of these contributions is given in Table 8. 
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Figure 4: Discriminant variable distributions for D* + — > D°7r + candidates selected in data 
(points with error bars) and in simulated hadronic Z decays (histogram) . 

Measured branching fractions, for right-sign events, are summarized in Table 9, where 
they are compared with similar results obtained in other experiments. The present analy- 
sis supersedes the previous DELPHI analysis [25] , in terms of statistical accuracy, as more 
information has been used to separate signal and background events, and systematics. 
All experimental results are compatible. 

There is no significant excess of events in wrong-sign combinations. Measurements of 
the D°7r~ and D + 7r + channels can be averaged, independently of the isospin of the irir 
system, as the same number of events is expected in the two channels giving: 

BR(b -> D°7i + 7i-rv e ) = BR(b -> D+Tr+Tr-r^) = (0.06 ± 0.08 ± 0.04) % (7) 

BR(6 -> D* + 7i + 7i-rv e ) = (0.00 ± 0.07 ± 0.07) % (8) 

The values corresponding to 90% C.L. upper limits are equal to, respectively, 0.18 % and 
0.13 %. 



Right sign 


External 


Tracking 


Lepton ID 


Fake lept. 


Discr. 


Cascade 


Total 


candidates 


BR,, Fib 


efficiency 


efficiency 


rate 


var. 


modelling 






0.028 


0.012 


0.016 


0.006 


0.047 


0.033 


0.067 


b^D+iv-rx 


0.021 


0.005 


0.006 


0.003 


0.007 


0.011 


0.026 


b^B*+n-£-X 


0.014 


0.007 


0.007 


0.003 


0.015 


0.013 


0.027 


Wrong sign 
















candidates 
















b -»• D u 7r-£-Jf 


0.028 


0.012 


0.016 


0.006 


0.047 


0.053 


0.079 


b D+Ti+rx 


0.021 


0.005 


0.006 


0.002 


0.007 


0.017 


0.029 


b^w+n+rx 


0.014 


0.007 


0.007 


0.004 


0.015 


0.067 


0.071 



Table 8: Contributions of systematic uncertainties (in 10 2 units) to the measured pro- 
duction rates ofb^ D^ttX£~Vi events. 
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b -> D u tt+£-X 


8.9 ±1.8 ±0.7 


11.6 ±2.4 ±1.1 


7.3 ± 1.8 ± 1.0 


b -»■ v> + TT-rx 


3.1 ± 1.0 ±0.3 


4.9 ±1.8 ±0.7 


3.0 ±0.7 ±0.5 


b D^TT^X 


4.0 ±0.9 ±0.3 


4.8 ±0.9 ±0.5 


4.7 ±0.8 ±0.6 


Wrong sign 








candidates 








b D u 7r"^-X 


2.6 ±1.3 ±0.8 


2.3 ±1.5 ±0.4 




6 -> D+7T+£-X 


-0.4 ± 1.0 ±0.3 


2.6 ±1.5 ±0.4 




6 -> D*+7r+£-X 


0.0 ±0.7 ±0.7 


0.6 ±0.7 ±0.2 





Table 9: Comparison between measured production rates (in 10 3 units) of 
b — > D^irXi'Ve events. 

3.4 Study of the D** hadronic mass distribution 

In order to study the mass distribution of right sign events, corresponding to the 
D** signal, the cut R > 0.25 has been applied on the discriminant variable to reduce 
the contribution from background events. This cut corresponds to a selection efficiency 
which varies between 67% and 85% for signal events, depending on the channel and on 
the data sample. 

A maximum likelihood fit has been performed using the A m = m(D(*V**) — m(D^) 
and the R variables, introducing the following components and parameters: 

• bo* = BR(B[j — > Dq + £~V £ ). In practice this component is considered to be a broad 
resonant mass distribution which can account for various possible states; 

• vtld*\ the mass of the Dq: it is kept fixed at 2.4 GeV/c 2 and a scan of its possible 
values between 2.3 and 2.5 GeV/c 2 has been made; 

• r D *: the total width of the D^; 

• b D - =BR(BQ ^D*+e-V e ); 

• m D *: the mass of the D^; 

• Yd*: the total width of the D^; 

• b Dl =BR(Bj->D+rF*); 

• 6 Da . =BR(B§_- D* 2 + £-V e ); 

• b^R — BR(B° — > T)Ti(rVe) corresponding to a possible non-resonant contribution; 

• Snr. the slope of an assumed mass distribution for the non-resonant component 
which is taken to be exponentially decreasing from threshold; 

• b n7T = BR(D** — > D7T7r) which is assumed to be independent of any particular D** 
state. 

Constraints from external measurements of the production rate and mass distribution 
for narrow states have been applied [33]: 

• b Dl = (0.64 ±0.11) %; 

• b D * = (0.28 ± 0.09) %. 

Removing these two constraints from the fit and imposing the values for the mass and 
width of the broad DJ state, we find b Dl = (0.33 ± 0.17) % and b D * = (0.37 ± 0.17) %, 
compatible with the values given above. 

The following decay channels have been considered: 



• D* -> Dtt; 
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• and Di — > D*7r and Dp; 

• the D?j can decay both into D7r and D*7r. The value 0.29 ± 0.07 [34] has been used 
for the decay probability into D*7r channels; 

• the possible non-resonant component is expected to decay into Dtc only as, for the 
D*7r channel, there is no contribution from the D* or B* poles which are expected 
to play the main role in non- resonant production [24]. 

As this analysis is based on the reconstruction of charged particles, decay channels with 
neutrals appear at a lower mass than the genuine decaying D**. Narrow states decaying 
into D*7r, where the D* decays into a neutral pion or photon and a reconstructed D, still 
appear as peaks in the D-zr mass distribution, but slightly broader and displaced from 
their nominal values. The expected accuracy of the mass reconstruction for completely 
exclusive decays and the effects induced by missing neutrals have been studied using 
simulated events. For completely reconstructed D** decays, the experimental resolution 
on A m is equal to 4 MeV/c 2 . For decays with a D* cascading into Dtc° or D7, and in 
which the neutral particle is lost, an additional smearing of 4 or 15 MeV/ c 2 is expected, 
respectively. For the D* + 7r~ exclusive decay channel, in which the D°, cascading from the 
D* + , decays into K~7r + (7r°), the mass resolution is poorer and values of 50 and 70 MeV/ c 2 
have been obtained from the simulation for 94-95 and 92-93 data samples respectively. For 
D** — > D7T7T decays, where only one pion is reconstructed, the Dtt mass distribution has 
been modelled using Gaussian distributions whose central value and standard deviation 
have been parametrized as a function of the central D** mass. Two parametrizations have 
been used for narrow and broad states, respectively. For broad states the variation of the 
standard deviation, as a function of the width of the D** state, has also been included. 
For narrow states with a 2.4 GeV/c 2 mass, the expected standard deviation of the D-zr 
mass distribution is 50 MeV/ c 2 whereas it is 90 MeV/ c 2 for a broad state having the 
same mass and a 200 MeV/ c 2 width. 

Events have been distributed in two-dimensional histograms of the discriminant vari- 
able R versus A m with a 10 MeV/ c 2 binning, to be sensitive to the presence of narrow 
components. For each channel (D°7r + (2) 7 , D + 7r~(l) and D* + 7r~(3)) and each data set (2), 
corresponding to twelve histograms in total, the observed number of events in each bin 
has been compared with the expectations assuming a Poisson distribution with average 
corresponding to the sum of the expected number of signal and background events. 

Two-dimensional distributions for the expected number of signal and background 
events have been constructed from the product of one-dimensional A m and R distri- 
butions. 

For background events, the A m distributions have been taken from simulation after 
verifying that their shape agreed with those measured with real data events selected in 
the sidebands (see Figure 5). A fit of these distributions, using the parametrization given 
by the following expression: 

Bckg{A rn ) = (A m - m n ) a exp (p + Pl A m + p 2 A 2 m + p 3 A 3 m + p 4 A 4 J; a = 0.5, (9) 

has been performed and the fitted distribution has been used to determine the expected 
average background in each bin. The R distribution corresponds to the one resulting from 
the fit of the different background components done in Section 3.3. The normalisation of 
background events has also been obtained from these fits. 

The A m distributions for signal events have been obtained by summing all components 
whose contribution is expected in the various channels. Breit-Wigner distributions have 
been used for each resonant state (i): 

7 The numbers in parentheses indicate the number of analysed samples corresponding to different D(°' + ) decay channels. 
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Figure 5: Comparison between A m distributions obtained for combinatorial background 
events selected in the sidebands for real data (stars) and simulation (diamonds) corre- 
sponding to the channel D°7r + in the 1994-1995 data sample. The two distributions have 
been normalized to the same number of entries. The curve fitted to simulated events has 
been superimposed. 



19 




possible displacement due to loss of neutrals in the D* decay, depending on the channel 
j of interest. The width, T^, receives contributions from both the natural width of the 
physical state Tj and from the experimental resolution. Expected variations of Tj, as a 
function of the mass (m) and of the angular momentum in the decay have been taken 
into account. Finally, iVj is the normalisation factor for the integral of BW(A m ) over the 
accessible mass range to be unity. 

For the non- resonant D-7T component, a normalized exponential distribution has been 
used: 



This distribution is maximum for low values of A m . Such a behaviour is expected when 
considering that the non-resonant Dir component is induced by the D* and B* poles, 
using chiral dynamics [24]. Using the same model, in the D*7r decay channel, the non- 
resonant contribution is expected to be small, as the D* and B* poles are not effective in 
this channel, and it has been neglected. 

The R distribution for signal events is taken from the simulation. 

3.5 Results on production characteristics of D** states 

Results of the fit are given in Table 10. The corresponding A m mass distributions are 
shown in Figure 6 and the fitted D** mass distribution, comprising all components, is 
given in Figure 7. 

Among the fitted parameters, three sets can be considered: 

• quantities corresponding to significant measurements in the present analysis: the 
production rate, mass and width of the broad state; 

• quantities introduced to parametrize components which are possibly contributing in 
the hadronic final state, but for which the present statistics do not allow a significant 
measurement. These components have been introduced as they can correspond to 
different measured D^7r mass distributions. For instance it is not equivalent to 
consider a measured Dn mass as originating from the non-resonant component or 
from a D7T7T decay. These parameters are the Dq production rate and width and the 
two quantities describing the possible non- resonant component; 

• quantities constrained by external measurements, in the fit, such as the production 
fractions of narrow D x and D2 states. 

In addition, two parameters (mo* and b nn ) have been kept fixed in the fit and then 
varied within a specified range to evaluate their contribution to systematic uncertainties. 
In the following some remarks on these results are made. 

3.5.1 Broad states Dg and 

Most of the D*7r component originates from the broad D\ state whose mass and width 
have been measured. These two last values can be compared with results obtained for 
this particle by the CLEO [35] and BELLE [36] collaborations, using B — > D^nrr decays: 



NR{A m ) 




(11) 



1 - exp (s NR (A. 



max 
•m 



m n )) 



20 





M(D*V)-M(D* + ) (GeV/c 2 ) 



Figure 6: A m distributions for right-sign events corresponding to the three D^n combi- 
nations. The fitted signal contributions are superimposed. 
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Figure 7: D** fitted mass distribution. The different components are shown within bands 
corresponding to the uncertainties of the fitted parameters. 



m D * = 2461^ ± 10 ± 32 MeV/c 2 ; T Dl = 290^ ± 26 ± 36 MeV/c 2 [35] 
m D * = 2427 ± 26 ± 20 ± 15 MeV/c 2 ; Yd* = 384±^> 7 ± 24 ± 70 MeV/c 2 [36] 

This is the first time that the parameters of this resonance have been measured in semilep- 
tonic B hadron decays. 

The present statistics is not sufficient to evaluate separately the contributions from the 
Dq (or other broad resonances), from the non- resonant component and from a possible 
D7T7T decay mode of D** states. 



3.5.2 The non-resonant component 

The fitted rate of the non-resonant component is compatible with zero and this is in 
agreement with expectations. The contribution from the D* pole to the non-resonant 
D7r component can be evaluated knowing the total branching fraction for the decay 

— > D* + £~Vi, the minimal value used for A m to select the events and the value of the 
D* width in the hypothesis that the mass is distributed according to a Breit-Wigner: 



D*po£e 

BR(B° -> D 7r + ri7 £ ) 



NR 



= ar ctan2 m(D07T+)lm ^ m ( D * + A x B R(B° - D*+r^) ~ 0.12% 

3 v 7r r D *+ j 

corresponding to a minimum D°7r + mass situated at 10 widths from the pole. For 
BR(B~ -> D+ 7r £ ^pJInr ' there is a natural cutoff corresponding to the sum of the 
D + and 7r~ masses and the value becomes 0.05%. 
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parameter 


Fitted value 
± stat. err. ± syst. 


b D{ 
m D * 


^l.z4 ± U.zo ± U.z ( J /o 
2445 ± 34 ± 10 MeV/c 2 
zo4 ± / 4 ± zo ivie v / c 


b D * 

bNR 
SNR 


(0.42 ±0.33 ±0.22)% 
260 ± 130 ± 130 MeV/c 2 
(0.23 ±0.35 ±0.44)% 
5±7 (GeV/c 2 )- 1 


b D *> 


(0.56 ± 0.10)% (constrained) 
(0.30 ±0.08)% (constrained) 



Table 10: Fitted values for parameters of the D** production characteristics. 



The contribution from the B* pole can be also included in a controlled way by consid- 
ering that the slow pion emission, by mesons containing a heavy quark, is governed by a 
universal parameter g whose value is fixed by the measured D* width [37]: 



r (D* - Dvr) = -£j- 2 pl (12) 



where f n is the pion decay constant and p n is the momentum in the D* rest frame. The 
interference term between the D* and B* pole contributions is also fixed and, using the 
values given above, this gives: 



BR(B° -> D°7i + £-V e ) 



D* + B*poles 

~ 0.07%. 

NR 



The main uncertainty, in the non-resonant contribution, comes from the interference 
with the Dq. Varying the branching fraction, BR(B^ — > D^i'Ve) between 0.1 and 0.8%, 
the non-resonant contribution changes between ±0.02% and —0.09% assuming that the 
interference has a negative sign below the Dq pole. The corresponding variation becomes 
±0.1% to ±0.2% for the other sign. 

3.5.3 Total D** production rate 

Summing all fitted components, the total rate for D** production amounts to: 

BR(B° -> D**rv e ) = (2.7 ± 0.7 ± 0.2)%. (13) 



3.5.4 Systematic uncertainties 

In addition to those considered in Table 8, systematic uncertainties have been evalu- 
ated for the following effects: 

• the b W7T branching fraction has been varied within the range (20 ± 15)%. The central 
value corresponds to an estimate obtained by comparing the phase-space integrals 
for the D*7r and Dp decay modes of a D** state. This variation is also in agreement 
with the value obtained in data for this parameter when it is left free to vary in the 
fit: (19 ±13)%; 
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• the central value of the Dq mass, m D *, which is kept fixed during the fit has been 
varied between 2.3 and 2.5 GeV/ c 2 ; 

• parameters used to define the shape of the combinatorial background distributions in 
A m (see Equation (9)) have been varied. The value of a has been changed between 
0.3 and 0.7 and the degree of the polynomial function was also taken as 3 or 5 ; 

• the contribution from B° semileptonic decays, which have been modelled similarly 
to the non-strange 6-mesons. Properties of narrow D** states published in [21] have 
been used. For broad states it has been assumed that their masses were displaced, 
relative to the corresponding non-strange states by the same amount as for narrow 
states. A m mass distributions for DMK + final states have also been evaluated by 
considering that the K + was reconstructed as a ir + . For excited charm states pro- 
duced in semileptonic decays of the A°, no simulation has been used, considering 
the lack of information on such decays. It has been assumed that the contribution 
from these final states was about 50% of the variation observed when B° decays are 
introduced. An uncertainty corresponding to twice this variation has been used to 
account for B° and A° contributions; 

• considering the recent result from the BELLE collaboration which has measured a 
larger width for D2 mesons [36] than quoted in [21], the effect of changing this value 
from 20 to 40 MeV/ c 2 has been evaluated. 

Individual contributions from these sources of systematics have been listed when eval- 
uating moments of the hadronic mass distribution. It has also been verified that possible 
additional systematic sources such as: 

- the relative branching fraction of D?; mesons into D*7r which has been varied accord- 
ing to the expected value: 0.29 ± 0.07 [34]; 

- the expected mass reconstruction accuracy which has been varied by 30% 

have negligible effects on hadronic mass moments. Uncertainties related to the control 
of the shape of the discriminant variable distribution have a small effect (4%) on the 
measured D** production rate and negligible contributions to hadronic mass moments. 

4 Moments of the hadronic mass distribution in b- 
hadron semileptonic decays 

Moments of the D** mass distribution can be evaluated from the results of the fit dis- 
cussed above. Results are given in Table 11. Statistical uncertainties have been obtained 
by propagating those on fitted parameters, using their full covariance matrix. System- 
atics are dominated by the uncertainty on the possible contribution from Dnn decays 
and could be reduced in future when experimental results on this decay channel become 
available. 

In determining the moments of the complete hadronic mass distribution in 6-hadron 
semileptonic decays, b — > D and D*£~z7^ channels have been included. 

For the first channel, values of branching fractions given in [21], have been used. As 
B° and B~ are expected to have the same partial decay width into the D£~z7^ channel, 
we get: 

BR(B° -> D + rv e ) = (2.06 ± 0.20)%. (14) 
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< TTlD** > 




V 771 £)* * ,> 






, 1U ^ 

777jj** 




CP„V /V 2s i 

^Vi6 V / C J 


^Vjre V / C J 


CPoV /V 2% i 4 

^Vi6 V / C J 


frvv /r 2 ^ 6 

yyj£ v / c j 


y\jt3 v / c j 


V / C J 


Veil 116 


2.483 


6.22 


40.1 


270.6 


1932 


14732 




±0.033 


±0.16 


±2.0 


±20.9 


±206 


±2039 




±0.030 


±0.14 


±1.6 


±14.4 


±122 


±1032 


m D , 
backg. param 


±0.008 


±0.04 


±0.3 


±2.5 


±18 


±36 


0.003 


0.02 


0.2 


2.6 


28 


291 


Bl, A" 


±0.010 


±0.04 


±0.5 


±4.3 


±38 


±334 


case, rate 


0.001 


0.01 


0.1 


0.4 


4 


25 


d± dist. 


0.002 


0.01 


0.2 


2.0 


22 


244 


T(D 2 *) = 40MeV/c 2 


-0.002 


-0.01 


-0.2 


-2.0 


-19 


-200 


Tot. syst. 


0.033 


0.15 


1.7 


15.7 


135 


1167 



Table 11: Measured moments of the D** mass distribution. When the sign of the variation 
is not given this is because the corresponding quoted systematic error originates from 
several sources corresponding to different signs which are given in the text. 



For the second channel, the value given in [38] at the time of the Winter 2003 confer- 
ences has been used: 

BR(B° -> D* + rv e ) = (5.27 ± 0.19)%. (15) 

The inclusive semileptonic branching fraction (BR(B[J — > ci~V()) has also been in- 
cluded as a constraint: 

T (W) 

BR(B° -> cTv t ) = -^BR(6 -> cTv t ) = (10.25 ± 0.30)%. (16) 

This last value has been obtained using the average values determined at LEP for the 
quantities entering into Equation (16), which are given in Tables 7 and 18. As the 
inclusive semileptonic branching fraction is the sum of the D, D* and D** contributions, 
the expected rate for D** production can be derived from these values: 

BR(B° -> D**£-V t ) = (2.9 ± 0.4)% (17) 

This is in agreement with the rate measured directly as given in Equation (13). 
Moments of the hadronic mass distribution have then been derived as: 

< rrijj >= Pd rrij-j + Pd* mJ^t, + pd** < m n D „ > (18) 

where pn = BR ^ B d~ >D ^ U£ \ value of pn.. has been obtained by imposing the con- 

straint 1 = po + Pd* + Pd** and including the measurement of the D** production rate 
as given in Equation (13). 

Results, following the notations in Equation (3), are given in Tables 12 and 13. System- 
atic uncertainties related to measurements of the branching fractions in Equations (14, 
16) are also given in Tables 12 and 13. 

Effects induced by the variation of the analysis efficiency versus the mass of D** states 
have been evaluated to correspond to an increase of 0.7, 1.5 and 2.5% for , Mif 1 and 
M% H respectively. Considering the present level of statistical and systematic uncertainty 
of actual measurements, these corrections have not been included in the quoted central 
values for hadronic moments. 

Error correlation matrices are given in Appendix A. 



25 





Mf 


Mf 


Mi 


Mf 


Mi 
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(GeV/c 2 ") 4 


(GeV/c 2 ) 6 


(GeV/c 2 ) 8 


(GeV/c 2 ) 10 


value 


0.647 


1.98 


7.4 


35.7 
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stat. uncert. 


±0.046 


±0.23 


±1.3 


±7.9 


±1080 


Ext. BR 


0.079 


0.22 


0.8 


4.1 


23.4 




±0.039 
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±0.6 
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±16.4 


m D . 


±0.015 


±0.04 


to.o 


T0.7 


T7.0 


backg. param 


0.007 


0.04 


0.2 


1.2 


8.0 


El, A° 


±0.007 


±0.03 


±0.2 


±0.9 


±5.2 


d± dist. 


0.005 


0.03 


0.2 


1.1 


7.4 


r(D%) =40MeV/c 2 


-0.004 


-0.02 


-0.1 


-0.9 


-5.8 


Tot. syst. 


0.090 


0.27 


1.1 


5.4 


32.3 



Table 12: Measured moments of the hadronic mass distribution, by reference to the spin 
averaged D-hadron mass, in b-hadron semileptonic decays. 
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±0.06 


±0.4 


±2.3 
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0.11 


0.7 


4.0 
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-0.10 


-0.5 
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Tot. syst. 


0.16 


0.32 


2.1 


11.7 



Table 13: Measured moments of the hadronic mass distribution, by reference to the aver- 
age mass squared, in b-hadron semileptonic decays. 

5 Moments of the lepton energy distribution in b- 
hadron semileptonic decays 

5.1 Inclusive reconstruction of 6-hadron semileptonic decays 

Selected events have been divided into two hemispheres using the thrust axis. The 
secondary hadronic system accompanying the lepton in the semileptonic decay has been 
reconstructed using an iterative procedure applied to the particles belonging to the same 
hemisphere as the tagged lepton. 

Charged particles, belonging to the hemisphere of the candidate lepton, with p > 0.7 
GeV/c, and with at least one associated hit in the Vertex Detector and positive impact 
parameter with respect to the primary vertex have been considered. These have been 
sorted in decreasing order of their probability of being B decay products based on their 
impact parameter significance and considering the particle that crosses the jet direction 
the furthest away from the primary vertex, or, when the crossing is not well-defined, the 
most energetic particle, to be the most likely B decay product. 

The charged particles have been iteratively tested for forming a secondary vertex. The 
procedure has been iterated while the following conditions have been fulfilled: invariant 
mass below 2.9 GeV/c 2 , distance from the primary vertex less than 3 cm but at least 
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2.8 times the uncertainty and on the positive side, and A% 2 < 3 after inclusion of each 
particle in the seed vertex fit. 

In those cases where no secondary vertex was found, single particles have been accepted 
when fulfilling one of the following criteria, in decreasing order of quality: a charged 
particle with p > 3 GeV/ c having a crossing point with the jet axis at least la downstream 
from the primary vertex and also downstream from the lepton candidate crossing point, 
but less than 15 cm from the primary vertex; a charged particle with p > 3 GeV/c with 
the largest impact parameter significance and positive impact parameter sign; the most 
energetic charged particle within 0.6 rad of the jet axis. 

The remaining charged particles in the hemisphere of the lepton with p > 0.5 GeV/c 
and at least one vertex detector hit not yet associated with the secondary vertex are 
then considered. Each of them is tested to belong to the vertex, and the one with the 
smallest x 2 contribution and giving a vertex mass closest to the D mass is included in 
the vertex, and all the remaining particles are then tested against this new vertex. This 
process is continued until the remaining particles have large contributions to the vertex 
X 2 and increase the mass difference between the vertex and the D mass. 

To improve the purity of the vertex, each of the particles associated with it is reconsid- 
ered, if there are more than two particles in the vertex, the xV n -d.f of the vertex exceeds 
2 and the mass of the vertex is more than the D mass. If removing a particle from the 
vertex improves the x 2 /n.d.f and the mass after the removal is closer to the D mass, the 
particle is permanently removed from the vertex. The average charged multiplicity of the 
vertices is 2.9. As the last stage, identified K®s within 0.8 rad and 7r°s within 1 rad of 
the lepton direction have been tested for association based on their energy, rapidity and 
contribution to the vertex mass. 

On average, 78% of the particles associated with the vertex were true D decay products, 
and 74% of the decay products were correctly associated with the vertex. The mass 
distribution of the reconstructed vertex is shown in Figure 8. 

The hadronic energy is obtained from the energy of the reconstructed secondary vertex, 
corrected with a linear function of the reconstructed vertex mass, when the reconstructed 
mass is below the D° meson mass. The direction of the hadronic system is taken from 
the momentum sum of the particles included in the vertex. 

For each decay, the energy of the B hadron has been estimated as the energy sum of 
the secondary hadronic system, the identified lepton and the neutrino energy. 

The energy obtained in this way was corrected by a function of the reconstructed 
hadronic energy, determined from simulation, with the maximum correction being 
±6 GeV. The neutrino energy was computed from the missing energy and momentum in 
both hemispheres corrected by a function of the missing mass in the event, determined 
from simulation. The resolution of the neutrino energy in B — > XiVi decays was estimated 
to be ±2.9 GeV. Neutrino energies larger than 1 GeV were required. 

The resulting resolution of the B energy was found to be 12% for 80% of all inclusive 
semileptonic B decays and 24% for the remaining decays. The resolution on the missing 
energy and on the reconstructed B hadron energy estimated with simulation for signal 
B — > Xciue events are shown in Figure 9. 

The direction of the momentum vector of the reconstructed XtPi system was adjusted 
by up to ±30 mrad with respect to the lepton direction using a function of the recon- 
structed B mass. Another estimate for the B hadron direction was obtained from the 
B line of flight reconstructed from the position of the vertex formed by the lepton with 
the identified charm charged decay products. An estimator, which combines these two 
independent measurements according to their expected resolutions as a function of the 
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vertex mass, GeV/c 2 

Figure 8: Mass distribution of the reconstructed charm vertex. The peaks at 0.14 an d 
0. 5 Ge V/c 2 correspond to vertices consisting of a single n or K , respectively. Points are 
the data, while the histogram is the simulation with the signal clear and the background 
hatched. 

reconstructed energies and the decay distance respectively, was defined. Resolutions of 
14(15) mrad have been achieved in (6) for 60% of all inclusive semileptonic B decays 
and 40 mrad for the remaining decays. 

The identified lepton was then boosted back to the reconstructed B rest frame and its 
energy E\ computed in this frame. This resulted in an average resolution of 170 MeV on 
E\ for 82% of all inclusive semileptonic B decays and 510 MeV for the remaining decays. 

5.2 Signal separation from background sources 

In the reduction of the b — > c — > £ and other backgrounds it is essential to avoid 
biases of the lepton energy spectrum. The separation was therefore performed using 
two discriminating variables, one based on the topology of the event and the other on 
charge correlations between the lepton and the other particles in the event, which are not 
sensitive to the lepton energy. 

The topological variable uses information on the lepton impact parameter with respect 
to the reconstructed secondary vertex, the topology of the tracks other than the lepton in 
the hemisphere, the number of particles not associated with the vertex in the hemisphere, 
the number of particles in the vertex and the y 2 of the vertex. 

The charge variable consists of a probability built from the correlation of the charge 
of the lepton and those of the reconstructed secondary vertex, of other vertices in the 
same and opposite hemispheres, of the jet charge of the opposite hemisphere and of the 
leading kaon candidate. The kaon candidate was identified based on kaon neural network 
output from the MACRIB package [39]. Two-dimensional distributions are shown in 
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Figure 9: The resolution on the missing energy (left), the fractional resolution on the 
reconstructed B hadron energy (centre) and resolution on the El energy (right) estimated 
with simulation for signal B — > X c £77£ events. 



Figure 10 for signal and background, respectively. The final separation variable (VSep) 
corresponds to the likelihood of an event to be a signal event, based on the location in 
the two-dimensional distribution of the topological and charge correlation variables. 

A further rejection of background is obtained by requiring a minimum value for the 
reconstructed B mass. In Figure 11 the distribution of the B mass is shown in data 
and simulation, for signal and background events. A B mass larger than 3.9 GeV/c 2 was 
required. A final sample of 14364 leptons was selected, with a purity in B — > X c £z/£ decays 
of 81%. 

Figure 12 shows the lepton spectrum, after VSep and B mass cuts, and the corre- 
sponding efficiencies as a function of E\. A lepton sample, depleted in signal by using 
an anti-cut on VSep is also shown, as a check of the shape of the simulated backgrounds 
with the data. 



5.3 Study of the lepton energy distribution 

The original lepton spectrum has been extracted from the reconstructed distribution 
by a spectrum re-weighting technique. This consisted of determining the resolution matrix 
relating the generated to the reconstructed spectrum for simulated signal events. Using 
this matrix, the coefficients of a re-weighting function for the generated spectrum have 
been fitted to minimize the y 2 between the resulting spectrum and that observed in the 
data. The efficiency correction has been taken into account at this stage. The procedure 
has been carefully tested on lepton spectra generated for different values of the and 
nl- parameters and smeared according to the resolution matrix. 

In order to increase the statistics in the signal description, the sample of simulated bb 
events (see Table 1) has also been used in the construction of the resolution matrix. 

A regularized unfolding method [40] has also been applied as a cross-check, but the 
re-weighting method has been preferred for its simplicity. 

5.4 Results on lepton spectra 

The resulting lepton spectrum is shown in Figure 13. 
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Figure 10: Two-dimensional distribution of the topological and charge correlation vari- 
ables for the signal B — > X C £T>£ (left) and the main background B — > X c — > Xiv (right). 



1600 




B mass, GeV/c 2 



Figure 11: Reconstructed B mass. Points are the data, histogram is the simulation with 
the signal clear and the background hatched. 
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Ej(GeV) E^GeV) 



Figure 12: The resulting El spectrum for samples enriched (left) and depleted (right) in 
B — > X c iVg decays using the separating variable. The data are shown as crosses while the 
simulation is shown as histograms. The different sources of background, shown shaded, 
are from bottom to top: b—>c^£,b—>-c^£, other lepton sources including decays 
in flight and converted photons, c — > t, misidentified hadrons. In the lower plots the 
efficiencies for the selection are shown, as a function of Eg. 




Figure 13: The background subtracted E\ spectrum (left) for the selected signal sample 
and the unfolded lepton energy spectrum (right). 
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The first, second and third moments, Mf, .\ /'/ and M' z , have been computed. In 
order to reduce the systematic uncertainties, the second and third moments have been 
computed with respect to the average value. 





Mt 




M'i 




(GeV) 


(GeV) 2 


(GeV) 3 


value 


1.3782 


0.1838 


-0.0301 


stat. uncert. 


±0.0073 


±0.0058 


±0.0015 


B species 


±0.0027 


=F0.0017 


T0.0005 


B -> D, D*, W*£u l 


±0.0010 


=F0.0005 


±0.0001 


B fragmentation 


±0.0027 


=F0.0020 


T0.0007 


b x u eu e 


±0.0008 


±0.0003 


to.oooi 


e.m. radiation 


±0.0035 


=F0.0001 


T0.0004 


Bkg modelling 


±0.0026 


to.ooii 


T0.0005 


B direction reconstruction 


±0.0027 


=F0.0018 


T0.0006 


B energy reconstruction 


±0.0027 


±0.0011 


T0.0003 


B mass cut 


±0.0051 


=F0.0031 


T0.0017 


Unfolding 


±0.0031 


=F0.0028 


±0.0029 


Tot. syst. 


±0.0092 


±0.0055 


±0.0036 



Table 14: Corrected lepton moments and sources of systematic uncertainties. 

The statistical correlation matrix for these three moments is given in Appendix B. 

In order to relate the measured moments to those computed for a B meson, some 
corrections need to be applied. 

Firstly the effect of radiation needs to be corrected for. This was done using correction 
factors computed separately for electrons and muons following ref. [41] and gives a shift of 
+7.0 MeV, -0.2 xl0~ 3 GeV 2 and -0.7x 10~ 3 GeV 3 on the first, second and third moments, 
respectively. Half of these shifts have been used as an estimation of the systematic 
uncertainty on this correction. Since e + e~ — > bb events at LEP result in the production 
of an admixture of 6-hadron species, a correction factor accounting for the bias due 
to the semileptonic decays of the heavier 6-hadrons was applied, using the fraction of 
B° and of 6-baryon left in the selected sample according to the simulation prediction. 
The uncertainty in the prediction is considered in the systematic uncertainty evaluation. 
Finally, the presence of b — > X M £z7^ decays results in a similar bias of the lepton spectrum 
toward higher energies, due to the larger phase space of this decay compared to b — > X c tp£. 
This was also corrected for. Results, after corrections, are given in Table 14. 

Several sources of systematic uncertainties have been investigated and the results are 
summarized in Table 14. The sources related to the modelling used in the simulation 
are: the fractions of the different D species in the decays, the different B species and 
the b fragmentation function. For the central value of the moments the same branching 
fractions for B — > D, D*, D** as in Section 4, the 6-hadron fractions from reference 
[21] and the results of reference [42] for the 6-hadron fragmentation distribution have 
been used, respectively. The variations quoted therein have been used for evaluating the 
systematic uncertainties reported in Table 14. 

The uncertainty related to the background modelling has been evaluated by changing 
the simulation prediction for the cascade decays within the uncertainties of the branching 
ratios given in [21] and changing the misidentification efficiency according to Sections 
2.2 and 2.3. It has also been checked by comparing the results obtained by using the 
background shape as extracted from the anti-tagged data with that using the prediction 
from the simulation. 
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The uncertainty due to the B reconstruction accounts for variations in the B recon- 
structed energy and direction. 

The uncertainty related to the unfolding procedure was evaluated by varying the 
reweighting function used in the fit and the binning. 

The stability of the results with respect to changes in the selection cut applied on the 
VSep variable and reconstructed B mass have been checked. Changes of the VSep cut 
inducing variations of the accepted statistics up to a factor of 1.5 and of the signal purity 
over a range from 76% to 89% have been considered and found to give results stable 
within their statistical uncertainty. The minimum value of the B mass required has been 
moved between 2 and 5 GeV, with a corresponding change of purity between 71% and 
86%. The maximum variation with respect to the central value has been used as an 
estimation of the systematic uncertainty in the agreement between data and simulation 
on the B mass reconstruction. 

Results obtained separately on the electron sample and the muon sample have been 
compared. The difference in the first moments amounts to 20±13±3 MeV, where the 
first uncertainty is statistical and the second is the uncorrelated systematic uncertainty 
on the background subtraction. The difference expected from simulation is 7 MeV. The 
differences on the second and third moments are fully compatible with the statistical 
uncertainty. 

6 Interpretation of the results 

A x 2 fit to the three leptonic (Mf, and M^, see Table 14) and three hadronic 
mass moments (M^ , M'^ 1 and M% H , see Tables 12 and 13) has been performed, using 
two theoretical frameworks. In the fit we also impose additional constraints derived from 
independent determinations. We follow the framework presented in [7] , updated with the 
new results discussed above, and using recent calculations given in [43] . Parallel fits have 
also been performed by other groups using several frameworks [44]. 

In the kinetic mass scheme, we fit the full set of parameters: m&(l GeV), m c (l GeV), 
p 2 and p z D . Expressions relating moments and these parameters can be found in [7]. We 
impose /ig=0.35±0.07 GeV 2 and plg-0.15i0.10 GeV 3 [43]. Two mass constraints have 
also been applied: m b (l GeV)=4.61±0.17 GeV/c 2 and m c (l GeV) = 1.14±0.10 GeV/c 2 as 
derived from the values quoted in [45], which were given using a different renormalization 
scheme. Results are obtained for a s (m&) = 0.22 ± 0.04 8 and are shown in Table 15. 
In order to study the effect of the bounds on m^c introduced, the fit has been repeated 
without these constraints. Results are shown in Table 16. Theoretical uncertainties have 
been evaluated following the procedure explained in [43]; namely 20% (30%) errors have 
been assumed for terms corresponding to 1 /m 2 (1/m 3 ) corrections and adding in quadra- 
ture variations corresponding to the uncertainty on a s {m h ). Corresponding theoretical 
uncertainties attached to m h and m c are due to those on a s {m h ). 

The central values of the heavy quark masses are in agreement with independent 
determinations [46,47]. The difference between the values of the two heavy quark masses, 
which are highly correlated, is m b (l GeV) -m c (l GeV)= 3.422 ± 0.034 ± 0.028 GeV/c 2 
(3.382 ± 0.051 ± 0.087 GeV/c 2 if no constraint on the quark masses is imposed). 

In the approach based on pole masses [48], the fit extracts A, Ai, A 2 , p\ and p 2 . We fix 
% = 0.0 GeV 3 and impose two constraints on Mb* — Mb and Mb* — Mb which effectively 
reduces by two the number of free parameters. The results are given in Table 17. 

8 For hadron moments a value of a s (nit,) = 0.3±0.1 has been used to account for missing terms in theoretical expressions. 
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Fit 


Fit 


Fit 


Syst. 


Syst. 




Parameter 


Values 


Uncertainty 


moments 


theory 




m b (l GeV) 


4.591 


± 0.062 


± 0.039 


± 0.005 


GeV/c 2 
GeV/c 2 
GeV 2 


m c (l GeV) 


1.170 


± 0.093 


± 0.055 


± 0.005 


GeV) 


0.399 


± 0.048 


± 0.034 


± 0.087 




0.053 


± 0.017 


± 0.011 


± 0.026 


GeV 3 



Table 15: Results of the fit in the m b (p), m c (p) and /j 2 (/j) formalism. 



Fit 

Parameter 


Fit Fit Syst. Syst. 
Values Uncertainty moments theory 


m b (l GeV) 
m c (l GeV) 
GeV) 

PI 


4.67 ± 0.11 ± 0.19 ± 0.03 GeV/c 2 
1.29 ± 0.17 ± 0.27 ± 0.04 GeV/c 2 
0.41 ± 0.05 ± 0.04 ± 0.08 GeV 2 
0.05 ± 0.02 ± 0.01 ± 0.03 GeV 3 



Table 16: Results of the fit in the m b (p), m c (p) and /i 2 (/i) formalism, without con- 
traints on m b (\ GeV) and m c (l GeV). Values of the fitted masses correspond to 
77^(77^)^ = 4.31 ± 0.20 GeV/ c 2 and M^frQ) 1 * 3 = 1.37 ± 0.24 GeV/ c 2 . 



Up to first order corrections in a s , parameters corresponding to non-perturbative QCD 
corrections, entering in the two approaches, are related: 

2 \ T ± + 3T 2 2 T 3 + 37^ 3 3 

p w = -Ai ; n G = 3A 2 + ; p D = pi; p LS = 3p 2 - (19) 

m b m b 

The parameter A enters in the expression relating heavy quark and heavy meson masses, 
which, for pseudo-scalar mesons reads: 

MB = mb + A + ^ b ~ + 4^ + °W b )> (20) 

where p 3 NL corresponds to a linear combination of 7i_4. Projections of the constraints 
from the six moments in the m b -p^ and m b -p 3 D planes are shown in Figure 14 and those 
in the A-Ai and A-pi planes in Figure 15. The xV n -d.f. of the fits is 0.4 and 0.2 in the 
two formulations. Since the contributions proportional to p\ s in the moment expressions 
are numerically suppressed, the fit is only marginally sensitive to its size and the result 
is determined by the constraint applied. By removing this, the fit would give p\ s = 
—0.4 ±0.4 GeV 3 . In contrast, the value of the leading 1/m 3 correction (parametrised by 
p%) can be determined with satisfactory accuracy and its range agrees with theoretical 
expectations [6] . These Figures illustrate the importance of the second hadronic moment 
to determine /j 2 (1 GeV) and of the second and third hadronic moments to extract p 3 D . 

In Tables 15 and 17, the first column of systematic uncertainties corresponds to sys- 
tematics on moments; correlated errors between the different moments have been included 
in the fit. The second column is due to systematics from theory. For the kinetic mass for- 
malism we propagate the uncertainty on a s and follow the suggestions of [43] to account 
for missing corrections. 
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r it 


r it 


r it 


byst. 


byst. 




Parameter 


Values 


Uncertainty 


moments 


theory 




A 


0.601 


± 0.065 


± 0.061 


± 0.05 


GeV 


Ai 


-0.252 


± 0.054 


± 0.018 


± 0.07 


GeV 2 


A 2 


0.117 


± 0.003 


± 0.002 


± 0.00 


GeV 2 


Pi 


0.032 


± 0.021 


± 0.010 


± 0.04 


GeV 3 


Pi 


0.085 


± 0.011 


± 0.017 


± 0.21 


GeV 3 



Table 17: Results of the fit in the A-Ai formalism. 



For the A-Ai formalism we include in the theory systematics the effect of % = (0.0 ± 
0.50) 3 , a s = 0.22±0.04 and we also estimate the effect of the missing corrections to third 
moments as explained in [7]. 




Figure 14: The projection of the constraints of the six measured moments on the 
mb(l GeV)-/i 2 (l GeV) (left) and m^l GeV)-p 3 D (right) planes. The bands correspond to 
the total measurement accuracy and are given by keeping all the other parameters at their 
central values. The ellipses represent the 1 o contours and include correlations between 
the parameters. 



6.1 Implications for \V c i,\ 

The value of |V C &| obtained from the total semileptonic decay width depends on the 
OPE parameters extracted above. We discuss now the implications of our results for 
|V c j,|, using the input parameters given in Table 18, which correspond to measurements 
obtained at LEP. 

The determination of |V C &| and the contributions of the various parameters in the 
kinetic mass scheme are described in [8]. This approach has been preferred to the pole 
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Figure 15: The projection of the constraints of the six measured moments on the A-Ai 
(left) and A-p± (right) planes. The bands correspond to the total measurement accuracy 
and are given by keeping all the other parameters at their central values. The ellipses 
represent the 1 a contours and include correlations between the parameters. 



Measurement 


Value 


Reference 


5-hadron lifetime 
BR(b -> Xfv) 
BR(6 -> XJru) 


1.573 ±0.007 ps 
(10.65 ±0.23) % 
(0.17 ±0.05) % 


HFAG Winter 2003 [38] 
LEPEWWG/2003-01 [49] 
PDG 2002 [21] 



Table 18: Input values, obtained at LEP, used for the determination of \V c b 



mass scheme as it does not rely on an expansion in l/m c and also because corrections 
contributing at order l/m\ have been fixed by experiment. Using the expression of \V c b\ 
quoted in [2] 9 , it gives: 

|Vcfrj = 0.0421 x (1 ± 0.014 meas . ± 0.014 fit ± 0.015 th .) , (21) 

where the first uncertainty reflects the accuracy on the semileptonic width determina- 
tion. This experimental uncertainty, corresponding to LEP-alone results, can be reduced 
to ±1% by including recent measurements obtained at B-factories which gave [38]: 
BR(6 — > X£~17f) = (10.73 ± 0.28)%, in agreement with the LEP measurement and hav- 
ing a similar accuracy. The second uncertainty, given in Equation (21), corresponds to 
uncertainties from the fit of the parameters, obtained in Table 15 and not including the 
theoretical uncertainties given in the last column. The third uncertainty has been taken 
from the result given in [2]. 

9 We have used the Equation (10) given in [2]. 
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7 Conclusions 

Production characteristics of D** mesons in 6-hadron semileptonic decays have been 
studied using exclusively reconstructed decay channels. 
The total production fraction has been measured to be: 



Decay final states are dominated by the D^*^ channel and no-evidence for a signal in 
channels with two pions has been obtained: 



BR(b -> D°7r + 7r-r^) = BR(6 -> D + n + n~rv e ) < 0.18 % at 90 % C.L. 



The dominant contributing channel is the broad whose mass and total width have 
been measured to be: 



Broad Dn final states favour a production which is maximum close to threshold, as 
is expected from non-resonant production, but the present statistics do not allow this 
feature to be firmly established. 

Moments of the hadronic mass distribution corresponding to D** states in 6-hadron 
semileptonic decays have been measured: 



< mo** 


> 


= 2.483 ± 0.033 ± 0.033 GeV/c 2 




> 


= 6.22 ± 0.16 ± 0.15 (GeV/c 2 ) 2 


4 

*\ TTt 


> 


= 40.1 ±2.0 ± 1.7 (GeV/c 2 ) 4 


*\ T7T/ j-)** 


> 


= 271 ±21 ± 16 (GeV/c 2 ) 6 


*\ 77% * 


> 


= (19.3 ± 2.1 ± 1.4) 10 2 (GeV/c 2 ) 




> 


= (14.7 ± 2.0 ± 1.2) 10 3 (GeV/c 2 ) 



Using these results and world averaged measurements for the exclusive 6-hadron 
semileptonic decay fractions into a D or a D* meson, various moments of the full hadronic 
mass distribution have been obtained in Section 4. 

Moments of the lepton energy spectrum in semileptonic B decays have also been mea- 
sured as: 



These results are interpreted in terms of constraints on the values of heavy quark 
masses, of the 6-quark kinetic energy and of the parameters contributing at order 1/m 3 
in theoretical expressions for the 6-hadron semileptonic partial decay width. The values 
obtained are: 



BR(B° -> D**rv e ) = (2.7 ± 0.7 ± 0.2)%. 



BR(6 -> D* + n + 7i-£-u e ) < 0.13 % at 90 % C.L. 



m Dl = 2445 ± 34 ± 10 MeV/c 2 
F D * = 234 ±74 ±25 MeV/c 2 . 



< E* > = 1.3782 ± 0.0073 ± 0.0092 GeV 

< (E* t - < El >) 2 > = 0.1838 ± 0.0058 ± 0.0055 GeV 2 

< (E* t - < El >f > = -0.0301 ± 0.0015 ± 0.0036 GeV 3 . 



m b (l GeV) 
m c (l GeV) 

GeV ) 




4.591 ± 0.062 ± 0.039 ± 0.005 GeV/c 2 
1.170 ± 0.093 ± 0.055 ± 0.005 GeV/c 2 
0.399 ± 0.048 ± 0.034 ± 0.087 GeV 2 
0.053 ± 0.017 ± 0.011 ± 0.026 GeV 3 , 
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and include corrections at order l/m\. 

Using these results, and inclusive measurements of the 6-hadron lifetime and semilep- 
tonic branching fraction obtained at LEP, an accurate determination of the value of the 
\V c b\ element has been obtained: 

|V c6 | = 0.0421 x (1 ± 0.014„ ± 0.014 fit ± 0.015 th .) . 

The first uncertainty becomes ±1% if measurements of B[J and B~ lifetime and semilep- 
tonic decay rates, obtained at the T(4S), are included. 

8 Acknowledgements 

We would like to thank P. Gambino and N. Uraltsev for helpful discussions during this 
work and for providing us with the necessary theoretical expressions. We had also useful 
discussions on the properties of the non-resonant D^tt component with A. Le Yaouanc. 
I. Bigi is warmly thanked for his guidance through the theoretical formalism. 

We are greatly indebted to our technical collaborators, to the members of the CERN- 
SL Division for the excellent performance of the LEP collider, and to the funding agencies 
for their support in building and operating the DELPHI detector. We acknowledge in 
particular the support of 

Austrian Federal Ministry of Education, Science and Culture, GZ 616.364/2-III/2a/98, 
FNRS-FWO, Flanders Institute to encourage scientific and technological research in the 
industry (IWT), Belgium, 

FINEP, CNPq, CAPES, FUJB and FAPERJ, Brazil, 
Czech Ministry of Industry and Trade, GA CR 202/99/1362, 
Commission of the European Communities (DG XII), 
Direction des Sciences de la Matiere, CEA, France, 

Bundesministerium fur Bildung, Wissenschaft, Forschung und Techno logie, Germany, 
General Secretariat for Research and Technology, Greece, 

National Science Foundation (NWO) and Foundation for Research on Matter (FOM), 
The Netherlands, 
Norwegian Research Council, 

State Committee for Scientific Research, Poland, SPUB-M/CERN/PO3/DZ296/2000, 

SPUB-M/CERN/PO3/DZ297/2000, 2P03B 104 19 and 2P03B 69 23(2002-2004) 

JNICT- Junta Nacional de Investigacao Cientifica e Tecnologica, Portugal, 

Vedecka grantova agentura MS SR, Slovakia, Nr. 95/5195/134, 

Ministry of Science and Technology of the Republic of Slovenia, 

CICYT, Spain, AEN99-0950 and AEN99-0761, 

The Swedish Natural Science Research Council, 

Particle Physics and Astronomy Research Council, UK, 

Department of Energy, USA, DE-FG02-01ER41155, 

EEC RTN contract HPRN-CT-00292-2002. 



; J ,8 



A Error matrices for hadronic mass moments 



In these matrices, which refer to moments Mf, Mf , Mf , Mf", Mf , Mf*, Mf*, 
and Mf* given in this order, the diagonal elements are the errors and non-diagonal 
elements correspond to correlation coefficients. 
Statistical error matrix: 



( 0.0455 



0.947 
0.232 



V 

Error matrix for systematics: 

/ 0.0896 0.981 
0.273 



0.865 


0.812 


0.921 





907 





725 


0.724 


0.698 


0.978 


0.951 


0.996 





994 





904 


0.899 


0.878 


1.29 


0.994 


0.992 





994 





973 


0.969 


0.956 




7.92 


0.974 





975 





990 


0.990 


0.983 






1080. 





998 





936 


0.933 


0.916 











176 





943 


0.936 


0.918 















736 


0.998 


0.993 
















4.53 


0.998 


















27.0 



\ 



V 



Total error matrix: 



0.949 


0.905 


0.863 





946 





655 


0.664 


0.540 


0.981 


0.945 


0.908 





990 





749 


0.747 


0.630 


1.09 


0.987 


0.967 





980 





848 


0.853 


0.758 




5.42 


0.995 





952 





910 


0.918 


0.844 






32.3 





920 





937 


0.950 


0.891 











161 





801 


0.789 


0.680 















318 


0.992 


0.974 
















2.06 


0.986 


















11.7 



/ 0.100 


0.945 


0.844 


0.759 


0.440 


0.872 


0.533 


0.544 


0.481 \ 


0.358 


0.966 


0.914 


0.664 


0.984 


0.763 


0.765 


0.712 






1.69 


0.987 


0.778 


0.987 


0.900 


0.903 


0.866 








9.59 


0.820 


0.957 


0.954 


0.959 


0.934 










1080. 


0.755 


0.870 


0.860 


0.851 












0.239 


0.853 


0.849 


0.804 














0.802 


0.997 


0.990 
















4.98 


0.996 


V 
















29.5 / 
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B Error matrices for lepton energy moments 

In these matrices, which refer to moments M[, \ //>' and M§ given in this order, the 
diagonal elements are the errors and non-diagonal elements correspond to correlation 
coefficients. 

Statistical error matrix: 

0.0073 -0.6041 -0.3435 \ 
0.0058 -0.5381 

0.0015 j 

Error matrix for systematics: 

0.0092 -0.7823 -0.2427 \ 
0.0055 0.0342 

0.0036 J 

Total error matrix: 

0.0118 -0.6942 -0.2578 \ 
0.0080 -0.1286 

0.0039 / 
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